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ABSTRACT 

We examine a new family of global analytic solutions for optically thick accretion disks, which 
includes the supercritical accretion regime. We found that the ratio of the advection cooling rate, 
Qadv, to the viscous heating rate, Q v is, i.e., / = Qadv/Qvis, can be represented by an analytical form 
dependent on the radius and the mass accretion rate. The new analytic solutions can be characterized 
by the photon-trapping radius, r trap , inside which the accretion time is less than the photon diffusion 
time in the vertical direction; the nature of the solutions changes significantly as this radius is crossed. 
Inside the trapping radius, / approaches / cx r n , which corresponds to the advection-dominated limit 
(/ ~ 1), whereas outside the trapping radius, the radial dependence of / changes to / cx r~ 2 , 
which corresponds to the radiative-cooling-dominated limit. The analytical formula for / derived here 
smoothly connects these two regimes. The set of new analytic solutions reproduces well the global disk 
structure obtained by numerical integration over a wide range of mass accretion rates, including the 
supercritical accretion regime. In particular, the effective temperature profiles for our new solutions 
are in good agreement with those obtained from numerical solutions. Therefore, the new solutions 
will provide a useful tool not only for evaluating the observational properties of accretion flows, but 
also for investigating the mass evolution of black holes in the presence of supercritical accretion flows. 
Subject headings: accretion: accretion disks, black holes — stars: X-rays 
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O 1. INTRODUCTION 

. 

Recent X-ray observations have revealed the existence of bright X-ray sources whose luminosities in some cases greatly 
"--■^ i exceed the Eddington luminosity of neutron stars, where the Eddington luminosity is given by Le = AttcGM/ k cs = 
43 1.25 x 1O 38 (M/M )erg s"\ where c, G, M, and k cs are the speed of light, the gravitational constant, the mass of 
the central object, and the electron scattering opacity, respectively. One of the best examples of such sources is the 
q ■ microquasar GRS1915+105 in our galaxy. This source emits mainly in the X-ray band, and with a luminosity of 
!— i ' Lx ~ 10 39 erg s _1 , in its high/very high state, close to the Eddington luminosity (Belloni et al. 1997, 2000). Narrow 
"o^ \ line Seyfert 1 galaxies (NLSls) are also candidates for near or super Eddington sources, with typical luminosities of up 
C\3 to 10 44_45 erg s _1 , despite their relatively small estimated black hole masses, M ~ 10 6 M Q (Bian & Zhao 2004; Collin 
L; ' & Kawaguchi 2004). 

. ,-h \ One of the more attractive explanations for such large luminosities is a supercritical accretion onto a black hole. The 
history of supercritical accretion theory is older than one might expect; for example, it is discussed by Abramowicz 
et al. (1978), Jaroszyhski et al. (1980), Paczynski & Wiita (1980), and others. However, most of the old models are 

5^ ! dynamically /thermally unstable (Begelman & Meyer 1982). One of the causes of an unstable solution is the lack of 
an advective cooling term in the energy equation. An advective flow for supercritical accretion transports the thermal 
energy into the black hole, so that the flow finally cools the disk and maintains it in a stable state. The first stable 
solution for supercritical accretion flow was proposed by Abramowicz et al. (1988), which is the so -called "slim disk" 
model (see in detail, Kato, Fukue, and Mineshige 1998, hereafter KFM98). Such a solution is seen as an optically thick 
flow with an accretion rate exceeding the critical rate, Merit = L^/c 2 — 1.3 x 10 17 (M/M Q )g s _1 — 2 x 10 _9 Af Q yr~ 1 . 
Several authors have developed this model and applied it to luminous objects (Szuszkiewicz et al. 1996, Wang et al. 
1999 for Active Galactic Nuclei (hereafter AGN); Mineshige et al. 2000 for NLSls; Watarai et al. 2000 for Galactic 
black hole candidates, hereafter GBHCs) where good agreement is found between predicted spectral properties and 
recent X-ray observations. Watarai et al. (2001) explicitly showed that a slim disk model can explain the thermal 
component of the spectrum observed in ultraluminous X-ray sources (hereafter ULXs). Ebisawa et al. (2003) performed 
direct spectral fitting of the slim disk model to observations of ULXs, and concluded that the slim disk model is the 
preferred model for the famous ULXs IC342 X-l. 

Simple models that can reproduce the exact numerical solutions are helpful to test the theoretical models observa- 
tionally. Self-similar solutions (hereafter SSS) for optically thin advection-dominated accretion flows (hereafter ADAF) 
were developed by Narayan and Yi (1994, 1995) and shown to be a good approximation to numerical solutions for a 
wide range of parameters (Narayan et al. 1997). Another class of self-similar ADAF solutions comprises the adiabatic 
inflow-outflow solutions (Blandford & Begelman 1999), which incorporate the effects of a jet or outflow. SSS of the 
convection-dominated accretion flow (CDAF) model have also been investigated by several authors who have compared 
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numerical results and SSS (Narayan et al. 2000; Quataert & Gruzinov 2000), and discussed their spectral properties 
(Ball et al. 2001). 

Self-similar solutions for the slim disk model were derived by Watarai and Fukuc (1999) and Wang and Zhou (1999). 
They explicitly derived effective temperature profiles of T c $ oc r" 1 / 2 , which agree well with those derived from the 
numerical results (Watarai et al. 2000). (In the case of the standard model by Shakura and Sunyaev (1973), the 
solution has already been expressed in analytical form, since the basic equations are not differential equations but 
algebraic equations.) However, it is not yet known whether other quantities, such as radial velocity, rotational velocity, 
and sound speed are reasonably well described by the SSS. One problem is that previous SSS differ from the numerical 
results shown in Figure 2 of Watarai and Mineshige (2001), and so we would like to find more useful solutions. One 
reason for this problem is that unlike the self-similar case in optically thin ADAFs, radiative energy loss plays an 
important role in optically thick solutions because the transition between the advection dominated regime and the 
radiative cooling dominated regime is very smooth. In contrast, the transition between an outer standard- type disk 
and an inner ADAF should be clear (see, e.g., Honma 1996). If such a solution can be found, it will be useful and 
sufficiently powerful to enable us to understand the flow properties, since it requires no numerical techniques. In this 
paper, we report new self-similar solutions that are a good approximation to numerical results over a wide range in 
radius, and we discuss the properties of the solutions. 

In the next section, we introduce the old SSS, following the same method as that adopted by Wang and Zhou 
(1999). In subsection E3 we derive the analytic form of / from the energy equation, and demonstrate the new 
analytic solutions by modifying the old solutions. In section we present profiles of physical quantities as a function 
of radius and compare the new solutions with numerical results. We then discuss the applications of the new solutions 
from the observational point of view in section The final section summarizes our conclusions. 

2. CANONICAL SELF-SIMILAR SOLUTIONS AND NEW ANALYTICAL SOLUTIONS 

In this section, we derive SSS based on the formulation of Wang and Zhou (1999) (see also Spruit et al. 1987; Narayan 
& Yi 1994). We assumed that the whole disk structure is optically thick. As in previous papers, we introduced the 
ratio of the advective cooling rate, Q~ dv to the viscous heating rate, Q^ is , f {— Qadv/Qvis)- 

In the optically thin ADAF, radiated photons can escape from the disk immediately, so that / keeps constant 
(/ ~ 1) throughout the whole of the disk. However, the ratio / should have a radial dependence in optically thick 
solutions. This is because / gradually changes as the radius decreases due to the radiative diffusion process. In the 
next subsection, we seek self-similar solutions without specifying the function of / as a first step. After that we will 
demonstrate the analytic form of / as a function of the radius and the accretion rate in subsection 12. 21 

2.1. Basic Equations and Canonical Self-Similar Solutions 

We solve the following hydrodynamic equations, (i.e. mass conservation, momentum equation in the radial direction, 

angular momentum conservation, hydrostatic balance), with the energy equation. Here we ignore the energy loss or 
the mass loss effects via an outflow. 

M = -27rrw r S, (1) 
dv r 1 dH , n0 n0 . ndlnil^ 

^ + E*= r ^-^-E-V £ ' (2) 

M{t-t in ) = -2^r 2 T rv , (3) 

(2N + 3)^ = H 2 n 2 K , (4) 

Qtis = *9adv + <9rad' ( 5 ) 

where S, £, £- ln , and H are the surface density , S = J podz, specific angular momentum (i — r 2 fl), angular momentum 
at inner edge of the disk, and the scale height. And we adopt the Keplarian angular frequency in a Newtonian potential, 
f^K = (GM/r 3 ) 1 / 2 . The r-<p component of the height-integrated viscous stress tensor was expressed by T rip , which is 
related to the total pressure by T rip — —an, where n is the height-integrated total pressure as defined by n = J pdz, 

and a is the viscosity parameter (Shakura and Sunyaev, 1973). We assumed a polytropic relationship, po = Kp^ +1 ^ N , 
in the vertical direction, where iV is the polytropic index. 

Viscosity generates thermal energy through the differential rotation of turbulently moving gas, and the viscous 
heating rate is expressed by 

Qvis = ~ rTrt f~^ ~ -T rv T n rt, (6) 

where Tq is the linear approximation form of angular velocity (Tq = — dlnfl/dlnr). Equation JJJ represents the "net" 
internal energy flux at each radius; i.e., advective cooling. 

M ds M n 
Qadv—^o— (7) 

Variables, To and sq represent the temperature and the entropy in the equatorial plane, and £ is a dimensionless 
quantity given by 

c _ 1 ( dlnpo r d\np \ 

^"fw 7hT7" r3 innTj' (8) 
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where I^ is the ratio of the generalized specific heats (e.g., Chandrasekhar 1967; Kato et al. 1998). Here £ is of order 
unity, and so we set £ = 1.5 for the remainder of this paper. 
We assumed the disk to be radiation-pressure -dominated, with opacity due to electron scattering only; 
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8cn 



2>npH k cs T,H 



(9) 



where k cs is the electron scattering opacity, k ~ k cs . 

To derive the characteristic radial dependence of the above equations, we transformed the equations 10, © 
as follows; 

n 



*2adv — fQvis 



(10) 



where B = (1 — £- m /£) is the boundary term. We also assumed f2 = r2 ^Kj where S7 is a constant. The momentum 
equations in the radial and azimuthal direction are 



n 



r 2 ^ 2 -^) = 0, 



(11) 

(12) 

where T v — d\nv r /d\nr, Tn = d\nll/d\nr, and the coefficients were determined from the SSS. We assumed Hq « 
BsjI and N = 3 for simplicity. Since the above equations can be solved analytically, we finally obtained the following 
solutions: 



Am Q n K B = 27raII, 



£1 = ^O^K; 

«7.17 x 10 3 
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1/2 



m 
100 
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:1 .0lxl0 7 / 1/2 (£) ~ 



10 



(13) 



(14) 



(15) 



(16) 



(17) 



The coefficients Tq, T v , and Tn can be determined self-consistently from the above equations (Spruit et al. 1987). 
In the limit / — ► 1, all of the solutions are consistent with Watarai and Fukuc (1999) and Wang and Zhou (1999). 
To summarize the normalized parameters, the mass accretion rate is m = M/M cr jt = M/(Le/c 2 ), the Schwarzschild 
radius is r g = 2GM/c 2 , and the black-hole mass is m = M/Mq, respectively. 
We assumed that the equation (fTfi)> is equal to the height-integrated radiation pressure, 



n = n r 



ad 



2/* g 



(18) 



where a is the radiation constant and /4(=128/315) is a numerical constant (Hoshi 1977). The radiative flux at each 
radius is 

F=7:Q; ad = -^ L =vT c V (19) 



3r 



Here Q rad is the radiative cooling rate in the optically thick medium (diffusion approximation). The optical depth 
r was defined by r = K es E/2 because the opacity is assumed to be electron-scattering-dominated. Hence, we can 



4 



New Analytical Formula for Slim Disks 



estimate the effective temperature distribution; 

1/4 



X,-, 



BTc 



/ V8 r l/4 fi V2 j 



K cs I 4 a y (27V + 3)£ 

1/4 / ~ \ V 



«2.52xf0 7 / 1/8 (^) f^j K (20) 

This temperature profile is characteristic of an optically thick ADAF. The slim disk behavior is expressed well by 
equation l|20l) . The radial temperature profile is flatter than that of the standard disk (changing from T c g oc r~ 3 / 4 
for standard disks to T c g oc r -1 / 2 for slim disks). It is interesting to note that the effective temperature does not 
depend on the accretion rate in this solution. However, the parameter / should be a function of the mass accretion 
rate because / ~ 1 corresponds to the high rh limit in the disk. That is, the effect of rh is implicitly included in /. 

2.2. Global Solutions including a Radiative- Cooling-Dominated Regime 

According to the numerical results of Watarai and Mineshige (2001), / profiles depend strongly on rh and radius. 
Strictly speaking, the energy equation used in the previous subsection is inconsistent because the adopted energy 
equation, Q~ dv = /Qvis' does not explicitly include the radiative cooling term, Q~ ad . The full energy equation should 
be Q~ dv = Qv is — Q ra d> an d so m the following analysis we start from this equation. 

First we included the effect of radiative cooling on / using the following form; 

Qvis Qadv + Qrad ^ 



where g is the ratio of radiative cooling to advective cooling. The explicit form of g is 

= Qrad _ W = 8 f-l/H. (22) 

■ Q~ dv Kcs HM£ y/(2N + 3)£Br n fio 
By inserting equation (1221 into equation l|21(l and l|17|l . we obtained a quadratic equation in /; 

f - (D 2 x 2 + 2)/ + 1 = 0, (23) 



where D and x are numerical coefficients (D = 8/y / {2N + 3)^BT n flo, for example, D w 2.18 for N = 3, and T n = 1.5), 
and there is a dependence on the radius and accretion rate {x = r /rh). This method is similar to an iteration procedure, 
and so we refer to it as a "semi-iterative method." 

It is clear from the discriminant of equation 12: ill that although two real roots exist for equation l|23|l , the range of 
/ is physically constrained to < / < 1. We can finally obtain an analytical solution for /; 



f(m, f) = f(x) = 0.5 D 2 x 2 + 2 - Dx\/D 2 x 2 + 4 . (24) 

In the limit x — > 0, f(x) tends toward unity, which corresponds to a high accretion rate disk. At the other extreme, 
i.e. x — > co, f(x) — > 0, and this case corresponds to a low accretion-rate disk. Note that in the above equation, 
/ is a function of the radius and mass accretion rates. This result is completely different from previous self-similar 
ADAF solutions. Figure ^ represents the radial distribution of /(m,f) for various accretion rates. The value of / 
goes to zero for smaller accretion rates and larger disk radii. For larger rh and smaller r, the / approaches to unity, 
i.e., becomes fully advection-dominated. These solutions are consistent with radiation- pressure-dominated disks. Note 
that the function / smoothly connects the advection-dominated regime and radiation-dominated regimes. When we 
replace / in equation (|14|l . 1)15(1 . 1(17|) with f(rh,f), the physical quantities of the advective-cooling-dominated regime 
are integrated with those of the outer radiative-cooling-dominated regimes. 

2.3. Photon- Trapping Radius 

For high accretion rates, some fraction of the photons that are generated in the disk cannot escape from the disk 
surface because they are trapped by the radial flow, so that the photons fall into the black hole along with the gas. It 
is natural that the efficiency of the conversion of gravitational potential energy to radiation energy decreases in this 
situation (Begelman 1978; Begelman & Meyer 1982). 

In the ID case, we can roughly estimate the radius of photon trapping, r trap , from a timescale argument. In order 
for photon trapping to occur, the accretion time must be shorter than the photon diffusion time, i acc < tdif, where t acc 
is the accretion time (~ r/v r ), t d a the diffusion time in the vertical direction (~ H/(c/t)), and r is an optical depth. 
We derived the trapping radius, r tra p, as follows; 



^*trap ~ 



o os |i/2 Mj (25) 
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Fig. 1. — The ratio of the advective cooling rate to the viscous heating rate. Solid lines represent the analytic expression of / in 
equation 1241 . By definition, the flow is advection-dominated above the dotted line (/ = 0.5). Input mass accretion rates are rh = 
10°, 10 ' 5 , 10\ 10 1 ' 5 , 10 2 , 10 2 ' 5 and 10 3 from bottom to top. 

so that the normalized trapping radius is, 



rtrap = — j{r tia . p ,m) 1 m. 
The previous equation is a biquadratic equation of ft rap , and so we can solve it analytically. 



^*trap — 
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\ 3D 2 



j51 
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loo y 



4 2| (3D 2 -4) 2 
(for D = 2.18). 



3D 2 



(26) 

(27) 
(28) 



Our analytic solutions roughly reproduce the numerical results not only for the inner advection-dominated regime 
(/ oc r° for r < r trap ) but also in the outer radiation-dominated regime (/ oc r~ 2 for r > r tra p) except in the 
inner/outer boundary regions (see figure |5J). 

Note that the trapping radius depends only on the accretion rate and not on the other disk parameters. That is, this 
radius is the transition radius between the radiative- cooling-dominated disk (standard disk) and the slim disk. The 
old SSS described in the previous subsection (i.e., the /=const. solutions) are valid only inside this radius. However, 
the derived SSS may not be reliable in the disk's inner region due to the condition used for its inner boundary. It is 
therefore necessary to be careful about the useful range of the solutions. 

2.4. Solutions in a Radiative-Cooling-Dominated Regime (r > rtrap) 
Outside the trapping radius, we simply adopted the standard solution of Shakura and Sunyaev (1973). Our solutions 



are basically the same as the standard ones so that the viscous energy input balances the radiative cooling, Q 
and the disk is radiation-pressure-dominated. 
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vis 



Q 



rad ' 



(29) 
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TABLE 1 

Summary of the radial dependence of the physical quantities 



Regime v r S T G ff H 

l(Q+ s ~Q a - dv ,p~ Prad ) ocr-Va ocr-i/2 « r" 1 ^ K r i 

Il(Q+ s ~Q- d ,p~p ra d) ocr" 5 / 2 ocr+ 3 / 2 oc r" 3 / 4 cx r° 

III (Q+ is ~Q~ ri ,p~p gaB ) ocr- 2 / 5 ocr- 3 / 5 oc r~ 3 / 4 ocr 21 / 



From the hydrostatic balance and angular momentum equations, we can rewrite the solutions as follows: 

A 2 AoO 3 



; —5.96 x 10 12 



2n \A{A 3 n% 



' 4.43 x 10 



)Q- < 32 > 



rn 

Toby 

. -I /4 / • \ 1/4 / x -3/4 



where Ai — BTq/(2tt),A 2 = 8c/n cs ,A 3 = (2N + 3)=9 (for A=3) are numerical constants. We set fl ~ B w 
1) To w 1.5 for simplicity. These solutions are all nearly the same as the Shakura and Sunyaev type solutions for the 
radiation-pressure dominated regime. Note again that these solutions can be adopted at r > r tra p. 

Here we summarize the radial dependence of physical quantities for each regime (see table y). In the following 
discussion, we will neglect the solutions for the gas-pressure dominated regime because such solutions have already 
been published in many papers. 

3. COMPARISON WITH NUMERICAL SOLUTIONS 

3.1. For Stellar Mass Black Holes 

Figure |21 shows a comparison of our new analytic solutions with numerical results. Numerical data are taken from the 
results of Watarai et al. (2000) and Watarai and Mineshige (2001). As is clear in figure[21 equation (|24|l reproduces the 
/ profiles of the numerical solutions relatively well except for the inner/outer boundary region. The trapping radius 
derived from equation (|24|l tends to be overestimated by a factor of two to three due to inner boundary conditions; 
however, using the new solutions, r tra p is roughly consistent with numerically obtained values for higher mass accretion 
rates. 

Note that the numerical solutions do not always correspond to an advection-dominatcd flow even if supercritical 
accretion takes place (m > 16, see also figure 2 in Watarai & Mineshige 2001). This shows that radiative cooling 
still works in the supercritical regime. Although this was demonstrated by Abramowicz et al. (1988), many authors 
do not recognize the fact that the slim disk is not always advection- dominated; i.e. Q~ ad is not always negligible. 
Therefore, when considering SSS in the regime of supercritical accretion, the radiative cooling term should be included. 
Figure shows the various physical quantities obtained from the new analytical solutions and numerical solutions for 
different accretion rates for a stellar mass black hole. The new solutions are indicated with thick solid lines, the 
numerical solutions with dashed lines, and the old SSS which correspond to / = const, are indicated by dotted lines 
for comparison. All of the analytic solutions are smoothly connected from the advection-cooling-dominated regime to 
the radiative-cooling-dominated regime. 

First, note that the velocity in the azimuthal direction roughly agrees with the new analytic solutions (top right 
panel in figure |2J). The deviation of v v between the analytic and numerical solutions is a factor of a few (/ ~ 0.1 
at most). The important point is that the differences from the numerical solution are a factor of 2-3, and the radial 
dependence of v v does not change to oc r -1 / 2 . This fact supports our initial assumptions that fl w (i.e., Sl ~ 1) 
and Tq 1.5 are valid over a large range of parameters. 

The radial- and sound- velocity, scale height, and effective temperature profiles obtained with the new analytic 
solutions are relatively consistent with those of the numerical solutions over a wide range of parameters. In particular, 
the effective temperature profiles are in excellent agreement with the numerical results for the whole disk structure for 
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Fig. 2. — Energy distributions of high rn flows. Solid lines represent the analytical form of / in equation 1241 . and dashed lines show 
the numerical calculations. The black-hole mass is fixed at m = 10. By definition, the flow is advection-dominated above the dotted line 
(/ = 0.5). 

higher mass accretion rates {rn > 10). This is because the effective temperature profile has a weak dependence on the 
inner-boundary-related quantities, i.e., T e g oc (BTn/S,) 1 ^ 8 (see equation I2(J|1 . 

Discrepancies between the analytic solutions and numerical solutions appear at the disk's outer region in all quanti- 
ties. This is because the outer region is in a gas-pressure-dominated state. Our analytical solutions assume a radiation 
pressure dominated state. Gas-pressure-dominated solutions have already been presented by Shakura and Sunyaev 
(1973) (see also Kato et al. 1998), and the boundary radius between the radiation-pressure dominated region and 
gas-pressure dominated region is located at 

f ab ~18(«m) 2 / 21 m 16 / 21 (34) 

/ a \ 2 / 21 /to\ 2 / 21 / rh \ 16/21 , . 

w601 (cu) Go) (loo J • (35) 

That is, if r > r a b, the disk is in the gas-pressure dominated state. Our new solutions are not useful in this regime. 

Numerical solutions for v r have steeper profiles around the inner boundary; however, the profile changes across the 
trapping radius. The standard radiation-pressure-dominated solution has a steep radial profile, oc r -5 / 2 , and so this 
profile agrees with the numerical solutions. Strictly speaking, the numerical solutions are expressed by the combination 
of the solutions oc r -1 / 2 and oc r~ 5 / 2 . 

The analytic solutions in the sound speed reproduce the solutions reasonably well. Here c s is derived from the 
relationship c s = (n/E) 1 / 2 . These results indicate that the thermodynamic quantities are expressed well by analytical 
solutions. These results can be explained by the fact that flow dynamics are influenced by radiative energy loss. 

The scale-height of the disk changes with radius, and our analytic solutions reproduce well the trend in disk structure. 
The disk geometry is moderately thick (slim) inside r trap ; outside r tra p, in contrast, the disk is geometrically thin. 
As for the surface density profiles (bottom left panel), the analytical solutions only reproduce the trend seen in the 
numerical solutions, but they differ slightly in normalization. The effects of the inner/outer boundary are significant 
for the surface density profile. 

As established in the classical model by Shakura and Sunyaev (1973), the a dependence of physical quantities related 
to the disk dynamics, density, surface density, and radial velocity is strong (p, E oc a -1 , w r oc «), whereas the quantities 
related to the temperature are not sensitive to the viscosity (dynamics), i.e., T e g, v v show no a dependence. (T c , c s , 
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Fig. 3. — Radial velocity (upper left), azimuthal velocity (upper right), sound speed (middle left), scale height (middle right), surface 
density (bottom left), and effective temperature (bottom right) profiles. The thick solid lines represent the effect of variation of accretion 
rates m=l, 10, 100, 1000 for our new analytic solutions, and the dashed lines are numerical solutions. The dotted lines show the variation 
of the parameter /=0.01, 0.1, 1.0 for previous SSSs. The dot-dashed lines represent the gas-pressure-dominant regime at each accretion 
rate. Velocities are normalized by light speed. 

and H weak have a weak a dependence.) We also show the a dependence of the radial and sound velocities in figure 
0] As mentioned above, the v r are significantly affected by the value of a, and the discrepancy between analytic and 
numerical solutions does not disappear, however, the sound velocities from the analytic solutions agree well with those 
obtained from numerical calculations even if a varies by three orders of magnitude. 

Note that these analytic solutions are not exact solutions , but represent the characteristic properties of the equations. 
However, the analytical solutions seem to be a good approximation for the temperature and sound velocity profiles in 
the disk's inner region, which is dominated by advective cooling. Therefore, the analytic solution is a useful tool for 
estimating the bolomctric luminosity in supercritical accretion disks. 

3.2. For Supermassive Black Holes 

As black hole mass increases, the size of the radiation-pressure-dominated region also increase. In figure |3] 
we demonstrate the supermassive black hole case. According to KFM98, the boundary radius, f a b, between 
the radiation- pressure dominated region and the gas- pressure dominated region is given by equation (|34|l : f a b ~ 
1800(a:/0.1) 2 / 21 (m/10 6 ) 2 / 21 (m/100) 16 / 21 for a supermassive black hole, i.e., an increase in black hole mass corre- 
sponds to an extension of the radiation-dominated region. Figure [S] shows that the analytic solutions for supermassive 
black holes are a better approximation than those for stellar mass black holes. However, it is necessary to consider 
carefully the effects of gravitational instability and self-gravity in supercritical accretion flows for supermassive black 
holes. We discuss these issues in section PTT1 



4. DISCUSSION 
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Fig. 4. — Radial and sound velocity for different values of the viscous parameter a = 0.001, 0.01, 0.1 from bottom to top. Mass 
accretion rate is set to be rh = 1000. Thick solid and thick dotted lines indicate the new analytic solutions for radial velocity, v r , and 
sound velocity c s , respectively. Dashed and small-dashed lines are numerical solutions. The dot-dashed lines around larger radii represent 
the gas-pressure-dominant regime. 

4.1. Gravitational Instability: Toomre's Q Parameters 

Gravitational instability in accretion disks is an important issue in understanding the mass supply mechanism from 
the outer regions of the disk. A great amount of study has been conducted in this area (e.g., Paczyhski 1978; Lin 
& Pringle 1987; Hure et al. 1998; Kawaguchi et al. 2004). In particular, Kawaguchi et al. (2004) investigated the 
effect of self-gravity in supercritical accretion flows and concluded that the origin of the thermal component of a NLS1 
observed in infrared band is the outer self-gravitating disk. In general, the gas density in a supercritical accretion disk 
increases so that the self-gravity effect may become important. The gravitational instability also depends on the value 
of a, i.e., if a is small (a ~ 0.01 — 0.001), the disk mass increases. 

To examine the gravitational stability of an accretion flows, we introduce Toomre's Q parameter (Toomre 1964), 

q = ~ 3.77 x KPom- 1 *- 1 / 8 ^ (36) 

where G is the gravitational constant. When Q is less than unity, the disk is gravitationally unstable. We can consider 
the Q — 1 radius; 

^-..»*io.(£n£)^(A)^(X). m 

This radius is much larger than the gas-pressure-dominant radius f a b, which means that the disk is stable to grav- 
itational perturbations in the radial direction. The Q parameters for 10 6 M Q are plotted in figure 6. The Q values 
for our new solutions are all larger than unity over a wide range in disk radius. Therefore, gravitational instability is 
not present the inner or middle disk regions for the range of disk parameters discussed here. Our result is supported 
by Kawaguchi et al. (2004) who also calculated Q values for NLSls, and found that Q is greater than unity for all 
disk radii. Note that for more massive black holes (M > 10 6 M Q ) or relatively low viscosity parameters (a < 0.1), the 
Q = 1 radius moves inward due to increasing gas density. (The Q value for stellar mass black holes is much larger 
than unity.) 



4.2. Total Luminosity 
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Fig. 5. — Same as Fig 2. but for the case where m = 10 . The viscosity parameter is set to be a = 0.1. 

Assuming spherical accretion, the maximum luminosity is determined by the balance between the gravitational force 
of gas and the radiation force. The mass accretion stops at the maximum luminosity, and this is the definition of the 
Eddington luminosity. In contrast, for disk accretion the gravitational force does not balance to the radiation force. In 
disk accretion, radiation escapes toward the vertical direction of the disk, so that the gravitational force cannot balance 
the radiation force. Eventually the mass accretion process leads to the production of large numbers of photons. The 
pressure gradient force (— is proportional to r~ 2 , which is the same radial dependence as gravity. The ratio of the 
gravity to the pressure gradient force is therefore always larger than unity, and hence the super-Eddington luminosities 
are possible. Numerical experiments indicate a maximum luminosity of ~ IOLe for rh = 10 3 — 10 4 (Abramowicz et al. 
1988; Watarai et al. 2000). 

We confirmed that our solutions are valid over a wide range of radii (see figure 7). In particular, the effective 
temperature from our analytic solution is in good agreement with numerical solutions. Accordingly, we can safely 
estimate the bolometric luminosity L at high accretion rates using our analytic solutions. A total luminosity fitting 
formula for supercritical accretion was first derived by Watarai et al. (2000). Fukue (2004) has expressed the analytical 
form of L using the self-similar formalism. Below we review again the self-similar determination of disk luminosity. 

The luminosity is mainly produced by two components; the inner slim disk Lslim and the outer standard disk £ssd : 

L = Lsjjm + issD 
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Fig. 6. — Toomre's Q for different accretion rates (m=l, 10, 100, and 1000 from bottom to top). In all cases, the black hole mass and 
the accretion rate are m = 10 6 and a = 0.1. Solid lines, dashed lines show the advection-dominated regime and the radiation-dominated 
regime, respectively. Dot-dashed lines represent the gas-pressure dominated solutions. 
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The bolometric luminosity normalized by the Eddington luminosity is 



L/L* 



4.0 J 'V 2 In 



4f t rap 



4f trap 



to > 20 
to < 20 



(38) 



(39) 



We divided the calculation into two regimes in order to avoid / < 0. The disk's inner edge was set to 3.0 r g for simplicity. 
The previous luminosity estimate derived by Wang and Zhou (1999) never exceeds the Eddington luminosity. The 
discrepancy between our results and theirs is due to a numerical coefficient and the second term of equation (|38|l . 
which is the contribution of a standard disk. The bolometric luminosity for various accretion rates is shown in figure 

In addition, we can calculate the bolometric luminosity using the new analytic solutions. The value of / is a function 
of the radius and the mass accretion rate, and so the /(f, to) should be included in the flux integration; 



L/Le = - 



BY, 



(2N + 3Kftg 



f(r,m) 1/2 f- 1 df 



(40) 



J 2.01 



f{f,rn) 1 ^ 2 f 1 df. 



(41) 



However, if we use equation l|41(l for small accretion rates (to < 20), the derived luminosities are significantly higher 
than the numerical results. This is because of the effect of the boundary condition for the inner region. As seen in 
figure |31 the difference between the numerical and analytical solutions increases for small m. Overall, equation l|39|l 
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Fig. 7. — Bolometric lu min osity as a function of mass accretion rate. The solid line represents the total luminosity using equation 1111 
for m > 20, and equation 1391 for small accretion rates (m < 20). The dashed and dotted lines are the numerical results and standard disk 
luminosity, L/Le = yj, respectively. 

gives a more accurate luminosity compared to equation (|41|) . If we consider the boundary condition more carefully, 
this discrepancy can be reduced. This will form part of future work. 

We note that the observed luminosity depends on the viewing angle. In the disk geometry, the emitting direction is 
mainly toward the vertical direction of the accretion disk. Thus, a super-Eddington luminosity will only be observed 
when viewing the disk face-on, whereas it will be measured to be sub-Eddington if viewed edge-on. It may be that the 
bolometric luminosity locally exceeds the Eddington luminosity in ID model. We should therefore take into account 
the effects of mass-loss to retain the self-consistency. 

The photon bubble mechanism is another promising model to produce super-Eddington luminosity (Arons 1992; 
Gammie 1998; Begelman 2002). The photon bubble instability leads to enhanced vertical transport of energy, so that 
the disk is efficiently cooled/or stabilized by the decreasing of thermal energy (Gammie 1998). However, numerical 
modeling of the photon bubble instability for high mass accretion rates has not yet been carried out (see Turner et al. 
2005 for subcritical rates). It is unclear whether the photon bubble mechanism can be effective in regions with large 
amounts of gas and high radiation densities. The effect of photon bubble processes in supercritical accretion flows is 
an important area of future investigation. 

4.3. Stability of Supercritical Accretion Flows 

The standard disk is always secular ly/thermally unstable in the radiation-pressure-dominant regime in classical 
studies (i.e., Lightman & Eardley 1974; Shibazaki & Hoshi 1975; Shakura & Sunyaev 1976). However, the stability 
analysis of this early work is only applicable to local stability; the early work assumed that the viscous heating energy 
is locally balanced by radiative cooling, and neglected the effect of advective energy transport, which is important in 
high mass accretion rate regimes. Therefore, these classical treatments donot simply apply to the present study. Global 
disk stability at high mass accretion rates has been studied by several authors. Linear stability analysis for high mass 
accretion flows indicates that the secular/thermal modes do not grow over a wide range in radius (Honma et al. 1991; 
Wallinder 1991a, 1991b). Moreover, the acoustic mode is also stabilized by advective motion (Fujimoto & Arai 1998). 
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Honma et al. (1991) also performed ID numerical simulation for relatively high mass accretion rates (M = lOAfEdd, 
where -Mgdd is the Eddington accretion rate, M^dd = 16Le/c 2 for pseudo-Newtonian potential), and confirmed that 
the middle (unstable) region of the disk is stabilized by the inner advection-dominated regime and outer gas-pressure 
dominated regime. Therefore, our new solutions are globally stable against secular/thermal instabilities. Recent 
2D radiation hydrodynamic (RHD) simulations by Ohsuga et al. (2005) showed that a quasi-steady disk structure 
appears for high accretion rates (e.g., rh= 100). Two-dimensional effects, e.g., convective motion, or outflows, may 
play an important role in preventing thermal instability. The photon bubble mechanism also stabilizes an unstable 
disk because it increase the radiative flux (Gammie 1998; Begelman 2002). Numerical experiments by Turner et al. 
(2006) have confirmed that the photon bubble mechanism produces an inhomogeneous density contrast, so that the 
disk can eventually radiate a large quantity of high energy photons from the disk equatorial plane effectively. It is 
difficult to full consider such a detailed scenario at this stage. In-depth comparisons with ID and 2D simulation results 
are required to confirm the stabilization mechanism of the thermal/secular instability. This will be the subject of 
future work. 

4.4. Is the Analytical Formula Useful for Observations? 

The greatest advantage of the analytical solutions presented here lies in their simplicity. However, analytical solutions 
may not be sufficiently accurate for detailed comparisons with observations because they do not incorporate the effects 
of boundary conditions. Nevertheless, we have seen that the global properties of the flow dynamics and spectra obtained 
from our solutions are roughly consistent with those obtained using numerical solutions. It is therefore worthwhile to 
discuss their observational implications. 

For radii less than the trapping radius, the temperature profile follows T e g oc j- -0 - 5 . and the spectral slope in this 
regime will follows vF v cx v . For radii larger than the trapping radius, the temperature profile of a standard disk, 
T e s oc r~ a75 , is expected, and the resultant spectral slope will be vF v cx is 4 / 3 . A spectrum observed to be flatter than 
the standard model would therefore be a useful tool for estimating the mass accretion rate. 

First, equation (|26|l shows that the trapping radius can be determined from the normalized accretion rate to. The 
temperature profile for high to disks changes at r — r trap ; therefore, the spectral slope must change at a characteristic 
frequency. If we can observe the turnover frequency in an accretion disk spectrum, we can estimate the accretion rate 
directly. The normalized photon-trapping radius is defined in this paper as r trap ~ to/2, and the effective temperature 
at rtrap is written as a function of /, m, and to. 

1.33 x 10 6 ^ UX" (%y 1/4 (f^r 172 for stellar-mass BH 
T eff (r trap ) » { ; f (i/8 _ 1/4 _ U o ( 42 ) 

7.50 x 10 4 K (gU (^) (fj%) for super-massive BH 



By adapting Wien's displacement law (i/ max = 5.88 x 10 10 T), and after considering the spectral hardening factor n 
(Shimura & Takahara 1995), we find the turnover frequency vto to be 

f 1.33 x 10 17 Hz (X) 1/8 (m)-V' ( )"V2 for steUar . mass BH 

Vto = "max W < ) ( 1/8 (43) 

I 7.49 x 10 15 Hz ( gfgj (Jjpr) 1 (f^) 1 (^) for super-massive BH 

For stellar mass black holes, the turnover frequency is located in the soft X-ray band (~ 0.54 keV). As the mass 
accretion rate increases, v io shifts to lower energies, so that it may be marginally possible to measure after considering 
soft X-ray absorption. According to previous work, the spectral hardening factor, k, may be larger than 1.7 in 
supercritical accretion flows (Kawaguchi 2003; Shimura & Manmoto 2003), in which case it will be possible to estimate 
the accretion rate robustly; however, it will also be necessary to carefully consider the effects of Comptonization. 

In supermassive black holes (M > 10 6 Mq) it may be possible to measure the turnover frequency. The position of 
the peak blackbody emission is located at ~ 10 14-15 Hz; which is within the soft UV/optical bands. Even if the black 
hole mass or accretion rate increases further, the turnover frequency vto still appears in the IR band; however, the big 
blue bump which is observed in the near-infrared band in quasar spectra cannot usually be observed directly due to 
the presence of strong hot dust emission (Chiang & Blaes 2003; Kishimoto et al. 2005). From a theoretical point of 
view, the effects of Comptonization, metal abundance, non-LTE, and a hot coronal component can all strongly affect 
to the big blue bump component in the quasar/AGN spectra (Hubeny et al. 2001). The determination of accretion 
disk parameters using the photon trapping radius may therefore be extremely difficult in practice. 

The origin of the high energy Wien spectrum is unknown due to the complexity of physical processes, for example, 
the relativistic Doppler effect, gravitational redshift, photon trapping, and Comptonization, which are expected to 
couple together. However, measurements of the turnover frequency could be a useful method to search for supercritical 
accretion candidates, or to obtain the information about the mass accretion rate. Using equation (|43|) . we can estimate 
the mass accretion rate, 

°ny' (wY 1 ' 2 (w^f ■ («> 

Unfortunately, the use of the photon-trapping radius to determine the disk parameters may be extremely difficult due 
to both theoretical and observational difficulties. 
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5. CONCLUSIONS 

1. We determined new analytic solutions to describe high luminosity accretion flows. Our analytic solutions are roughly 
consistent with the global numerical solutions for the same assumptions (except at the inner or outer boundary). In 
particular, we explicitly showed that the ratio of the advective cooling rate to the viscous heating rate / is a function 
of the radius and mass accretion rate. 

2. Thermodynamic related quantities obtained using the analytic solutions, in particular, the temperature distribution 
gives a reasonable approximation to the profile of the whole disk. This demonstrates that our solutions are useful 
for interpreting observations, although complications may arise due to emission and scattering in the disk corona and 
possible Compton scattering. 

3. Using observations of the flux ratios of particular objects, we can discriminate between a standard disk and a 
slim disk. In other words, this is an indirect method of measuring the accretion rate. For this measurement, the 
key physical quantity is again the photon-trapping radius. Hence, the estimation of the trapping radius is of greatest 
theoretical importance. Our model is a simple ID model, and so it necessary to calculate the exact trapping radius 
using more than 2D or 3D radiation hydrodynamics. 

In the present study, we did not consider the effect of outflows. In high mass accreting flows, the gravitational 
force in the vertical direction of the disk cannot balance the flux force, and in some cases, the flux locally exceeds 
the Eddington limit. Hence, the vertical structure of supercritical accretion flows will also be an important topic for 
future studies. 

Two-dimensional RHD simulations are technically complex and Time consuming tasks. Our new analytic solutions 
for optically thick disks will therefore be useful in developing a new understanding of luminous black hole candidates. 
We expect to use the future observational data to investigate the supercritical accretion scenario in more detail. 

We would like to thank R. Narayan, S. Mineshige, J. Fukue, H. Kamaya, K. Ohsuga, and A.K. Inoue for helpful 
comments and discussions. The author also would like to thank C. James for checking the manuscript, and thanks the 
Yukawa Institute for Theoretical Physics at Kyoto University, where this work was initiated during the YITP-W-01-17 
on "Black Holes, Gravitational Lens, and Gamma- Ray Bursts". This work was supported in part by the Grants-in 
Aid of the Ministry of Education, Science, Sports, and Culture of Japan (16004706, KW). 
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